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Abstract 

The Bethe-Salpeter equation is related to a generalized quantum-mechanical Hamil- 
tonian. Instability of the presumed vacuum, indicated by a tachyon, is related to a 
negative energy eigenstate of this Hamiltonian. The variational method shows that 
an arbitrarily weak long-range attraction leads to chiral symmetry breaking, except 
in the scale-invariant case when the instability occurs at a critical value of the cou- 
pling. In the case of short-range attraction, an upper bound for the critical coupling 
is obtained. 
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1 Introduction 



There are many important contexts in which chiral or other symmetries are dynamically 
broken due to a condensation of (say) fermions. It is accepted that in QCD with massless 
quarks, the colour-singlet flavour-singlet quark bilinear has a non-zero vacuum expecta- 
tion value resulting in a spontaneous breaking of chiral symmetry. Moreover, this symmetry 
is expected to be restored at high temperature. Other condensates of fermion composites 
are also suggested as the baryon chemical potential is varied. Many exotic condensations 
are suggested in technicolour theories. Condensation of relativistic massless fermions is also 
important in the context of effective field theories for many condensed matter systems. 

The most popular way of handling chiral symmetry breaking is to phenomenologically 
incorporate it through an effective Lagrangian. Almost the only a priori calculation is 
through a nonperturbative solution of the Schwinger- Dyson equation with some truncations. 
This was pioneered long ago by Nambu and Jona-Lasinio pQ. This was in turn motivated 
by the gap equation in the theory of superconductivity. 

The BCS theory of superconductivity [2] is founded on the Cooper pair instability. The 
instability of the normal vacuum can be demonstrated for an arbitrarily small attraction 
(mediated by the phonon) between the non-relativistic electrons. It may be noted that it 
is easier to establish this instability than to calculate the consequent condensation. 

We would like to have a criterion for instability of the presumed vacuum for the case of 
relativistic massless fermions in quantum field theory. The analysis of a fermion-antifermion 
pair in this case is complicated by various issues in the theory of Bethe-Salpeter equation 
(BSE). We consider Euclidean BSE at zero energy-momentum, and relate the equation to 
a generalized quantum-mechanical Hamiltonian. Instability indicated by a tachyon pole 
in the Bethe-Salpeter amplitude is then related to a negative energy eigenstate in the 
quantum-mechanical system. 

It is to be noted that this approach does not involve solving the gap equation. Our 
calculation is performed with massless fermions, and a tachyon appears at some critical 
value of the coupling. This indicates instability of the presumed, chiral symmetric vacuum. 

We now list some other works which consider the issues addressed in this paper. The 
instability of the normal state in the BCS theory is demonstrated in Ref. (3j using ladder 
diagrams for electron-electron scattering. The criterion for instability of the chiral sym- 
metric phase in QED corresponding to the existence of tachyonic bound states of fermion- 
antifermion pair was considered in Ref. jl] (see also the review in Ref. [5]). The gap equation 
was related to a quantum mechanical Hamiltonian to find the critical coupling for the onset 
of chiral symmetry breaking in Ref. [B] for quenched QED and in Ref. [7] for gauge theo- 
ries with extra dimensions. The ladder approximation for the scalar vertex in QED3 was 
reduced to a quantum mechanical problem in Ref. [Sj. Ref. [9] addresses critical couplings 
in a variety of situations in condensed matter systems. 

The paper is organized as follows. In Sec. [21 we consider the case of massless QED 3 
and explain our approach for establishing the instability of the presumed vacuum. In Sec. 
[31 we introduce the use of the variational method for locating the critical coupling, and 
illustrate the method for massless QED3. In Sec. HI the variational method is applied to 
treat the generic situation with massless exchange, including both the scale-invariant and 
the scale-noninvariant cases. The case of short range attraction is addressed in Sec. In 
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Figure 1: The Bethe-Salpeter equation at vanishing total energy- momentum for the 
fermion-antifermion pair (the box represents the 2PI scattering kernel). 




Figure 2: The scalar channel leads to a drastic simplification in the spinor algebra. 
Sec. El we present our conclusions. 



2 The case of massless QED3 

To establish a proper setting, we first consider the case of massless QED in three (space- 
time) dimensions which has many interesting features. Vacuum polarization due to iV 
massless fermion flavours in a 1/N expansion changes the infrared (IR) behaviour of the 
photon from inversely quadratic to inversely linear in momentum [IDl [TTJ. This feature 
survives to all orders [TT], |T2]. Moreover in a specific non-local gauge [T3"j [Ti] . there are 
no logarithmic IR contributions to the fermion self-energy and the vertex correction to all 
orders [12]. Thus there is no anomalous dimension for the fermion in this gauge, and the 
full vertex can be replaced by the bare vertex in the IR. 

Therefore we consider the BSE in the ladder approximation: this is given by Fig. [T] 
with the simplest kernel in the form of a single photon exchange. Let us now consider the 
following special situation [15] . We work with the Euclidean theory. The photon propagator 
is replaced with its IR behaviour, which suffices because any instability is expected to be 
driven by the attraction at small momentum transfer. We consider fermion-antifermion 
scattering at vanishing total energy- momentum and in the scalar channel (see Fig. [2]). This 
provides a drastic simplification in the spinor algebra, and we are rid of all the spinor and 
Lorentz indices. We then obtain a simple integral equation: 

| q — q ' | J q" z \ q — q "\ 
Here the coupling A equals 8/37r 3 iV to the lowest order in 1/N [15J. We may convert this 



integral equation into a differential equation by applying V^: 

( " Vj - ^jmn = ^\q- q') . (2) 

Thus I(q, q') is the Green function for the Hamiltonian of a three-dimensional Schrodinger 
equation (in variables q) with an attractive inverse square potential 

V{q) = -Ait\/q 2 . (3) 

Therefore 

n e " 

where Xn{q) are the eigenf unctions of the Hamiltonian with eigenvalues e n . This is a well- 
studied problem in quantum mechanics with unusual properties [16] . When the attraction 
is smaller than a critical value, A < A c = 1/167T, the Hamiltonian has a continuous spectrum 
starting at e = 0. For A > A c , the spectrum is unbounded from below. 

As we are considering the Euclidean BSE with a vanishing total incoming momentum, 
e n < in Eq. (J3J) corresponds to a tachyonic bound state of the fermion-antifermion pair. 
This signals instability of the system against chiral symmetry breaking [15] . A study of 
the gap equation by other authors has led to a dynamical mass formation at precisely 
this critical value of the coupling [13J. The important thing to note is that the approach 
simplifies the problem in so much as relating it to a negative energy eigenstate for a quantum 
mechanical Hamiltonian. This simplification is utilized below. 

In this model, there is a sudden change to arbitrarily negative energy values at a A c . 
This is an artifact of this particular model: exact scale-invariance (in the IR) leads to a 
continuum of negative energy eigenvalues, once there is one such eigenvalue. 

Note also that the inverse square potential is coming from the massless fermion, and 
not from the forces of attraction. The photon propagator gives rise to the 'kinetic energy' 
operator (— being the operator for which the propagator l/\q—q'\ is the Green function). 



3 The variational method 

In the case discussed in Sec. [2} we were fortunate in having an exactly solvable quantum- 
mechanical Hamiltonian to locate the critical value of the coupling. This is not the situation 
in the other cases to be discussed below. We apply the variational method to locate the 
critical value of the coupling at which a tachyon is formed. We first illustrate the technique 
for the case of massless QED3. 

We want an estimate of the ground state energy. For this we should use a trial wave 
function which is large at q = so that the potential energy of Eq. has an expectation 
value (V) (which is necessarily negative definite) as large in magnitude as possible. In 
addition, the wave function should not fall off too rapidly in order that the the expectation 
value (T) of the kinetic energy (necessarily positive definite) is not too large. Let us 
therefore use a trial state 

iP(q) = q- s e -^ (5) 
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with s > and [i > 0, using which we can make exact calculations. Normalizability of this 
wave function requires 3 — 2s > 0. Requiring that (V) is finite gives a stronger condition 
1 — 2s > 0. In that case, (T) is also finite. We have 

||^|| 2 = 4vr(2 / i) 2s - 3 r(3-2 S ), (6) 



(V>|V|V> = -167r 2 A(2/i) 2s - 1 r(l - 2s) . (7) 
As the wave-function has no angular dependence, V 2 ip is effectively (l/q)d 2 (qip)/dq 2 and 

(^|(-V 2 )|^) = 2 7 r(2/i) 2 - 1 (l - S )r(l - 2s) . (8) 



Adding Eqs. ([Z]) and (jHJ) and dividing by 



(E) 



2s 



2 , we obtain 
8ttA 



(9) 



Since 1 — 2s > 0, and therefore, 1 — s > 1/2, the quantity within the square brackets 
cannot be made negative for A < 1/16-7T. Whereas for A > 1/16-7T, one can increase /i to get 
arbitrarily negative value for (E). We thus get the correct value for the critical coupling. 
The change of the minimum possible value of (E) from zero to — oo as A crosses A c is a 
consequence of the scale invariance in the problem. 



4 Generic situation with massless exchange 

We now apply our technique to a generic situation with massless exchange (Fig. [1]). We 
work with the renormalized theory. Presuming the chiral symmetric vacuum, we consider 
the IR behaviour after all IR divergences in this theory have been taken into account. Let 
the leading behaviour of the 2PI fermion-antifermion scattering kernel in the scalar channel 
for small momentum transfer \q — q'\ — ■> be K,/\q — q'\ a - (As earlier, we consider the case 
of zero total incoming momentum.) We also allow an anomalous dimension for the massless 
fermions, and take the fermion propagator to be j\j q^ +1 in the IR. Thus the BSE (keeping 
only the leading IR behaviour) in d space-time dimensions is 

m n = w^ff + A / ddq, '^w^r I{(^> ' , " ,) • (10) 

Here A = ac/ (27r) d , and we have factored out one factor of k from each term of the equation. 
We now obtain the formal operator for which the Green function is l/\q — q'\ a . Consider 
the Fourier transform 

\q-q'\° C[,a) J(27TY e W-a- ^ 

The required operator is therefore (l/c(d, a))(V 2 p d ~ a ^ 2 where = —id/dq^. We may 
regard q as the 'coordinate' variables and p as the conjugate 'momentum' variables, with 
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the corresponding operators Q and V (with % = 1). We can now formally relate the integral 
equation (jTUl) to a generalized quantum-mechanical Hamiltonian: 

^=(^ 2 ) ( ^ Q)/2 -Ac(rf,a)^, (12) 



HI(q,q') = c(d,a)5W(q-q'). 



(13) 



As in the previous case, we want to know when there is a negative eigenvalue for this 
Hamiltonian. It suffices to choose the same trial wave-function as given in Eq. ([5]). The 
calculation of (T) is easier when d is an odd integer (because the angular integration is 
simple). Therefore we present the calculations only for d = 3. To calculate (ip\T\ip), we use 
the wave- function ip{p) in the 'momentum' variables, since T is diagonal in the momentum 
basis: 



d 3 qe ip ^(q) 



271 T(2 s)( 1 1 

ip — ip) 2 ~ s (fi + ip) 2 ~ 



We then obtain 

(tP\(V 2 f- a)/2 \ 



d 3 p 



2„3-a 



2[r(2 - s)] 2 fi 2s - a 
a a 



x 



B[2--,--s) -2cos((2--)7r)£?(4-a,a-2 S 



(14) 



(15) 



(16) 



where we put p = [i tan 9 when doing the radial integration and used the integrals given in 
Eqs. (EED and ([29]). Also, 

(i/>\V\if>) = -47rAc(3,«)(2 / u) 2s+2/3 - 3 r(3-2s-2/?). (17) 

Adding Eq. (|T6|) and Eq. (jTTJ) and dividing by H^H 2 (as given by Eq. ([6])), we obtain (E) 
for the general case in three dimensions. 

We first consider the scale-invariant case. As a first example, it can be checked that for 
the QED 3 case a = 1 and (3 = 1, the above general result for (E) correctly reproduces Eq. 
([9]). (The formula given in Eq. (131]) is useful for this purpose. Note also that c(3, 1) = 47r.) 
Next consider the case a = 2, (3 = 1/2 in three dimensions. This case is also scale- invariant. 
Now we have to use c(3, 2) = 2n 2 . We find 



(E) 

Here s < 1 from finiteness of ( 
for < s < 1 is 2/ ^pK for s = 
after Eq. ([9]). 



2vr 2 /i 



1 r(3-a) 
2vr 5 /2r(5/2-s) 



A 



;is) 



0|T|^) and (V>IW)- Th e smallest value of T(3-s)/r(5/2-s) 
1. This gives A c = 1/n 3 using similar argument as presented 



6 



Let us now consider the general scale-invariant case d — a = 2(3 in d space-time dimen- 
sions, so that T and V scale in the same way. By scaling q — > q/ fi, we immediately arrive 
at the form 

(E)=^f 1 (s)[l-\/f 2 (s)]. (19) 

Here f\(s) and f 2 (s) are positive definite functions of s in the allowed range of values of 
s. (This is because (T) is positive definite and (V) is negative definite.) Following the 
earlier arguments, we conclude that the minimum value of f 2 {s) gives the critical coupling 
A c . Again (E) can be made arbitrarily negative by taking fi to be arbitrarily large when 
A > A c . 

We next consider the scale-noninvariant case d — a ^ 2(3. Consider first, for example, 
the case a = 1, (3 = 1/2 in three dimensions. In this case the general formula gives 

For any given A, one can choose /i < 47rA(l — 2s)/(l — s) with any s in the allowed range 
0<s<l/2to have a negative value for (E). Thus there is instability for any coupling. 
(Note that on completing the square in Eq. (I20I) . one finds that the minimum possible value 
of (E) is — 4ir 2 X 2 , which is finite as expected in the absence of scale invariance.) 
For the general scale-noninvariant case in d dimensions, 

(E) = 9l (a)fi a - \g 2 (s)v b (21) 

where gi(s) and g%(s) are positive definite functions of s, and a ^ b. For a > b, one can 
choose any s in the allowed range and [i < (\g 2 (s) / gi(s)) 1 ^ a ~ b ' to have a negative (E) for 
any given A. Similarly for b > a, one only needs to chose \i large enough. Thus, chiral 
symmetry breaking takes place for an arbitrarily weak attraction in the scale-noninvariant 
case involving massless exchange. It may be noted that this is the situation for QCD in 
four space-time dimensions. 



5 The case of short-range attraction 

Finally, let us address the consequences of a short-range attraction. To present the main 
features, we limit ourselves to a massless fermion of canonical dimension and an attractive 
force provided by the free propagator of the mediator. Now the BSE is 

J (£ ?') = F 4rT 2 + A / d V4 r 4rT 2 J (^ ?') • ( 22 ) 
\a — a r + m z J a z \a — ar + m z 



Since 



1 2 7 r 2 / 7 ^e-(-')^, (23) 



\q — q'\ 2 + m 2 J (2ir) 3 p 

the formal operator for which this is the Green function is given by 

T ~ v^e mv/ ^ (24) 
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Note that the spectrum blows up exponentially. 

We cannot take the earlier trial wave-function because (ip\T\if)) does not exist: ip{p) 
falls off only as power of p for large p. T is not an operator on the Hilbert space of such 
wave-functions. We need to choose a different trial wave-function precisely because the 
spectrum of T blows up so fast. Indeed we cannot choose wave-functions which are singular 
near q = 0, as we did earlier to gain as much potential energy as possible. Any power law 
behaviour of i/j(q) as q —>■ translates to a power law behaviour of ip(p) as p — > oo, and 
then (ifj\T\tp) will not be finite. This is the reason that it is so difficult to have negative 
energy eigenvalues now. We choose the simple trial wave-function 

V»(p) = (25) 

p 

which corresponds to 

^ ) = w + /.') ' (26) 

so that we can make explicit computations. Thus (ip\T\ip) is finite if \x > m/2 and we 
obtain ll^ll 2 = l/(4vr 2 /i), (^\T\ip) = 1/(2tt 2 (2// -m) 2 ) and (i;\V\ip) = -A/(2ii 3 ). Therefore 

/i 3 - 47r 2 A(/i - m/2) 2 
(E) = 2^ -m/2) 2 • (27) 

Now /(//) = /i 3 — 47r 2 A(;U — m/2) 2 is positive both at /x = m/2 and for large /i. It is negative 
for a range of values of /i in between, provided A is sufficiently large; then /(//) possesses 
a minimum. This minimum value of f(fi) is precisely zero when A = A c . Therefore setting 
both /(//) and df/dfj, equal to zero for A = A c , we obtain A c = 27m/327r 2 . Since other 
choices of i(>(q) may give a lower value for A c , this is actually an upper bound for the critical 
coupling. 



6 Conclusion 

In this paper, we presented a criterion for instability of the vacuum against chiral symmetry 
breaking. The instability, as indicated by a tachyon pole in the Bethe-Salpeter amplitude, 
was related to a negative energy eigenstate of a quantum-mechanical Hamiltonian. Our 
approach was illustrated for chiral symmetry breaking in massless QED 3 , and the variational 
method was then applied to address the critical coupling for generic situations. 
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Appendix 

Some of the integrals and formulas used in this paper are listed below [TH 



o 



" 12 dx sin" 1 " 1 x cos"" 1 x = -B I — , - I (28) 



n-i 1 D fm n\ 
2 V 2 2/ 



/•7r/2 TTi7T 

/ fix sin™ - - 1 x cos™ -1 x cos(m + n)x = cos B(m, n) (29) 

jo 2 



_B(m, n) = r(m)r(n)/r(m + n) (30) 



2 2m-l / j\ 

r(2m) = — — r(m)r m + ~ (31) 



7T V ' \ 2 
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